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Abstract. We consider factorizations of the stationary and non-stationary Schrodinger 
equation in M™ which are based on appropriate Dirac operators. These factor- 
izations lead to a Miura transform which is an analogue of the classical one- 
dimensional Miura transform but also closely related to the Riccati equation. In 
fact, the Miura transform is a nonlinear Dirac equation. We give an iterative 
procedure which is based on fix-point principles to solve this nonlinear Dirac 
equation. The relationship to nonlinear Schrodinger equations like the Gross- 
Pitaevskii equation are highlighted. 
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1. Introduction 

Also today nonlinear equations and operators are a challenge. But it is also 
well-known that transformations do a good job in understanding and solving non- 
linear equations. One famous transform is the Miura transform ([27], 1968) which 
in one spatial dimension represents a connection between Korteweg-de Vries equa- 
tions. The Miura transform can also be considered as a special case of the Riccati 
equation and can be used to factorize the (stationary) Schrodinger equation. In 
higher dimensions Clifford algebras and analysis give the possibility to generalize 
the factorizations into higher spatial dimensions and even to space-time domains 
(t, x) G M+ x G C K+ x R n . 

Factorization of second order differential operators have been studied several times. 
Factorizations of the Helmholtz operator in case of a real wave number had been 
studied by K. Giirlebeck p] and K - Giirlebeck/W. SproBig in their books [12]. [T5|. 
Z. Xu [33], [34j, F. Brackx and N. van Acker [6] and together with R. Delanghe 
and F. Sommen [2j considered generalized Dirac operators of type D + k, with 
the Dirac operator D and A; is a complex number. E. Obolashvili [2E!,[29] and 
later on Huang Liede [14| investigated the case of purly vectorial k. Fundamental 
solutions for the Dirac operator D + a, where a is a paravector, are described in pj. 
Factorizations of the Helmholtz equation and the related Dirac operators D + a, 
where a is a complex quaternion are investigated by V. V. Kravchenko and M. V. 
Shapiro [20], [21]. In the context of Maxwell's equations operators generalized 
Dirac operators considered in [24] by A. Mcintosh and M. Mitrea. 
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After this introduction we give a short overview of the Clifford operator calculus in 
general, starting with Clifford algebras and then describing Clifford analysis and 
operator calculus as they needed for applications to equations of mathematical 
physics independently of a specific Clifford algebra. 

In the next section examples are presented, where the main emphasize is laid on 
the Dirac operator and associated monogenic functions. The last example are 
parabolic Dirac operators which allows to factorize the heat equation. 
Section 3 is concerned with the non-stationary Schrodinger operator and appro- 
priate factorizations by parabolic Dirac operators. 

In Section 4 we investigate the factorization of the stationary and non-stationary 
Schrodinger equation which will lead to the Miura transform. 
With the Miura transform as a nonlinear Dirac equation and a procedure based 
on fix-point principles to solve this nonlinear equation is section 5 concerned. 
The last section deals with the Gross-Pitaevskii equation and it's relation to a 
nonlinear Schrodinger equation and the Miura transform. 

2. Clifford Operator Calculus 

The connection of analysis and algebra has created an deeper understanding of 
physics and natural sciences. On way to algebra to create analysis can be done by 
the use of Clifford algebras. 

2.1. Clifford algebras. Let e 1} . . . , e n be an orthonormal basis of K n . The Clifford 
algebra G? n ,o is the free algebra generated modulo 

2 I 1 2 

x = —\x\ eo, 

where e is the identity of the Clifford algebra and the basis elements fulfil the 
anti-commuting relations: 

where <5y denotes the Kronecker symbol. Let N = {1, 2, . . . , n) and for each 
non-empty subset A of N set 

e A = e Al . . . e Ak , A = {A u A k }, 1 < A x < A 2 < . . . < A k < n. 

By convention e$ = e the identity of the Clifford algebra. 

There are three operations on a Clifford algebra. The principal automorphism 

e' A = (-l) |A| e A , AC N, 
where \A\ is the cardinality of A, the principal anti- automorphism 

e A = (-l)^-% A , 
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while the conjugation is the decomposition of the principal automorphism and 
anti-automorphism 

e A = (ei)' = (e> A y = {-1)MM+% A . 

It is easily seen that 

G'i — G"i ■ % — X ^ • . . j TJj •. 3.H.d G<i G j — Gj G^ • 

We will identify the Euclidean space W l with /\* Ci Qtn the space of all vectors of 
G?o,n, i.e. any element iGl" can be identified with 

n 

X ^ ^ Xi&i. 
i=X 

In this way it is easily seen that each non-zero vector has inverse given by ™. 
More on Clifford algebras and Clifford analysis can be found in the pioneering 
book |4j, the treatment of elliptic boundary value problems is considered in [12j. 
The function theoretical basis and the application to more complicated partial 
differential equations like the Navier-Stokes equation is contained in [TJ|. A (com- 
prehensive) function theory for the Dirac operator is the book [5j. A good de- 
scription of Clifford analysis and its connection to classical as well global analysis 
can be found in [9]. 

2.1.1. Special cases. 

1. Complex numbers C ~ d ,i '■ Identifying 1 with e and the complex unit i with 
e\. We obtain that the complex numbers C are equivalent to the Clifford algebra 

(Xo,i. 

2. Quaternions HI ~ G? ,2 : Let us denote the basic elements of the quaternions by 
1, i, j, k. Then 

ij = k = —ji, jk = i = —kj, ki = j = —ik, 

and we identify 1 = e , % — e ly j — e 2 and hence k = eie 2 . Thus the quaternions EI 
can be identified with G? ,2- 

Remark 1. Generally Clifford algebras Cl VA are defined by the generating 
elements % — 1, . . . , e„, n — p + q, where 

= -1, i = 1, ■ ■ p, e} = +1, i =p + 1, . . . ,p + q = n, 

and eiej + = —25ijeo. 

But we restrict ourself to the cases Ci 0jn and Cl n $ because it is possible to 
factorize all typical types of second order differential operators by using these 
Clifford algebras. 
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2.2. Function spaces. 

Let G C K™ be a bounded domain with sufficient smooth boundary T = dG, whose 
complement contains a non-empty open set. Then, functions / denned in G with 
vales in (X , n are considered. These functions may be written as 

/o) = X^ AeA ' x e G ' 

A 

where = • ■ ■ = e^.../^, 1 < hi < . . . hk < n, and e0 = e . Properties such 
as continuity, differentiability, integrability and so on, which are ascribed to / have 
to be possessed by all components /a (a?)- In this way the usual Banach spaces of 
these functions are denoted by C a , L p , W£. The norm in the Banach-spaces B is 
given by 




In particular, we introduce in L 2 (G) of real-valued function f A the G? ,n-valued 
inner product 

(/, g) ■= / f{x) ■ g(x) dx, 
Jg 

where f{x) = J2a Ia^a- The inner product is related to the L 2 -norm by 

l|/|li» = 8C</, /)=^IM 2 (X), 

A 

where scg = g e denotes the scalar-part of g = J2a9a- 
2.3. The operators D, T, F. 

In this section we describe in a general way the basic operators which will become 
specified by the use of a specific Clifford algebra. In general, let G be a domain 
with smooth boundary dG = T and denote by n(y) the outward pointed normal 
at y G T. There are two ways of consideration, the classical one which assumes the 
domain G to have a Ljapunov boundary T, i.e. T is of class C 1 ' with arbitrary 
a > [l2j, [T2j. But these strong assumption on the boundary can be weakened 
to the consideration of strongly Lipschitz domains [22], [23j. 

We consider a differential operator D which is a " square root" of a second order 
differential operator and therefore possesses a fundamental solution e(x), i.e. 

De(x) = 5(x) in G. 

Because the multiplication in Clifford algebras is not commutative we have to 
distinguish between left-monogenic functions, i.e. functions / such that Df = 
and right-monogenic functions g, which fulfil gD = 0. Then the convolution 
integral operator 

Tu = —e-k u = — I e(x — y)u{y) dy 
Jg 
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is a right but not a left inverse of the operator D and there exists an operator 

(Fu)(x) = J e(x - y)n(y)u(y) dT y , x E G, 

such that 

DT — I and TD + F = I. 

The operator F appears by applying the formula of partial integration. By taking 
limits we obtain 

_ lim Fg = ±g + Sg, 

inside G/outside G^r 

with a Hilbert-type operator 

(Su)(x) = 2 J e(x - y)n(y)u(y) dT y , x e T. 

2.4. Decompositions. Another general principal is that 

L 2 (G) = ker£>©im£>*. (1) 

It is a challenge to prove that the (orthogonal) complement is imD*, this set is 
closed itself. Now, it is easy to see that this decomposition defines (orthogonal) 
projections P and Q. 

Remark 2. The just mentioned decomposition are orthogonal in L 2 and can 
be extended to direct sums in If, 1 < p < oo, and in Sobolev spaces W£, p > 1. 

The connection between the operator F and the Hilbert-type operator S give 
raise to another decomposition which is based on the (orthogonal) projections 
P = \{I + S) and Q = \{I — S). This is a Hardy- type decomposition 

L 2 (T) = imP©iniQ. (2) 

An interesting connection is given between ((H) and (j2j) in the way that / G ker D 
if and only if tr / 6 im P. This connection is given by the Plemelj-Sokhitzkij's for- 
mulas which tell us that tr / e im P if and only tr / has a monogenic continuation 
/ into G, i.e. Df = 0. 

2.5. Examples. 

2.5.1. Dirac operator. Here, D denotes the Dirac operator 

3=1 J 

which factorizes the Laplacian: 
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where A n denotes the Laplacian in M n , and we obtain the fundamental solution 
of the Dirac operator by applying the Dirac operator to the fundamental solution 
of the Laplacian H. We get the generalized Cauchy kernel 

-1 x 
e[x) = — — — , 

where u n denotes the surface area of the unit sphere in R" . The first conclusion is 

Corollary 1 (Borel-Pompeiu formula). Let G be a bounded strongly Lipschitz 
domain ofR n . If f G WHG), 1 < p < oo, we have 



Ff + TDf 



f, inG, 
0, inR n \G. 



Similar to the case of one-dimensional complex variables we can define Hardy 
spaces [25j. Let us denote by G + = G and (7_ = M. n \G. 

Definition 1. For a function F : G± — > Ci ,n o,nd 1 < p < oo we set 

\\F\\ HP := sup { [ \F{x ± S)\>dT x Y , 
<5>o Ur J 

and define the Hardy spaces of monogenic functions 

W(G±) := {/ is left monogenic in G±; ||/||kp < oo}. 

There is also a connection to the boundary values of functions from the Hardy 
spaces H P (G±). 

Corollary 2. Let 1 < p < oo. For a left monogenic function H in G, the 
following are equivalent 

(1) H e H P (G); 

(2) There exists h e L P (Y) such that H is the (right) Cauchy integral 
extension of h, i.e. H(x) = Fh(x) for x in G; 

(3) H has a non-tangential boundary limit H + (x) at almost any point 
x G T, i.e. there exists 

H + (x) := lim H(y) for a.e. x G T, 

and H is the (right) Cauchy integral extension of its boundary trace. 
In addition the norms \ \ ■ \ \ HP and \ \H + \\ LP are equivalent, i.e. 

\\H\\jq> = 1 1 H 1 1 . 

The most important singular convolution operator is the singular Cauchy inte- 
gral operator 5 on T, defined for almost all x G T by 

Sf(x) = 2 lim / e(x - y)n(y)f(y) dV y . 

£ ^°J{yer:\x-y\>e} 
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Corollary 3 (Plemelj-Sokhotzkij's formula). Let G be a bounded strongly Lip- 
schitz domain o/M n with boundary dG = T. If f e L P (T), 1 < p < oo, t/ien 

(Pf)(x) = lim / e(x-y)n(y)f(y)dT y + -f(x) 

=l(I + S)f(x). 

and hence 

{Qf){x) = \{I-S)f{x). 

Moreover, we have 

Sf = 2Pf(x)-f(x) and Sf = -2Qf + f. 

The equality of the integrals is shown in [22]. It follows from Calderon-Zygmund 
theory. We have 

Corollary 4. The space I/ P (r), 1 < p < oo, posses the decomposition 

L P (T) =imPnL p (r)©imQnL p (r), 

where the decomposition is orthogonal in case p = 2 and direct in case p ^ 2. 
Moreover, the spaces imPfi L P (T) and imQ n L P (T) can fee identified with the 
boundary values of the Hardy spaces W(G±). 

The Teodorescu transform is given by 

Tf(x) = - [ j^-f(y)dy 

u n J G f - y\ 

and from Borel-Pompeiu's formula we immediately conclude that 

trTf G imQ. 

2.5.2. Cauchy-Riemann operator. Now, D denotes the Cauchy-Riemann opera- 
tor 

Df = K + y e .K and Df = ^--Ye^- 
dx ^—f 3 dXj dxn ^ 3 dXj ' 

3 =l J j=l J 

Also the Cauchy-Riemann operator factorizes the Laplacian. We have 

DD — A n+ i, 

where A n+ i denotes the Laplacian in R n+1 , and we obtain the fundamental solution 
of the generalized Cauchy kernel 
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where u> n+ i denotes the surface area of the unit sphere in W n+1 . 
The Cauchy-Riemann operator can be transform into a Dirac operator by multi- 
plication with e n+ \ : 



dxn ^ 3 dxj n+1 dxn n+1 ^ 3 n+l dx.j dxn n+1 ^ 3 dXj 
j=l J j=l J j=l J 

On the other hand the Dirac operator can be transformed into the Cauchy- 
Riemann operator: 

^ 63 dx, { 6n) ~ dx n + ^ { ^ dx, ~ dx n + ^ j dx, ' 

2.5.3. d-operator. The simplest case occurs when the underlying Clifford algebra 
is just equivalent to the complex numbers, i.e. (3? ,i ~ C. In this case we have the 
complex Cauchy-Riemann operators 

« 9 , • 9 a a 9 ■ 9 

9 * = TT +l TT and °' = a 

ox oy ox oy 

which can be identified with 

^ d d - d d 

D = ^ + ei TT and D= ^~ ei ^~ 
ox oy ox oy 

where 5 is identified with D and we have the properties 

- - d 2 d 2 

dd = 3d = DD = DD = A 



dx 2 dy 2 

The fundamental solution of the <9-operator is the Cauchy kernel 

11 11. 1 z 

e(z) = = 1 ~ -. — tttCi 

K J 2m z 2ixz 2n\z\ 2 

and we have the generalized Cauchy formula (Cauchy-Green formula, Borel for- 
mula, ...) as an anolog of the Borel-Pompeiu formula: 

2vr i J r C ~ z tt J g ( - z 

2.5.4. Generalized Dirac and Cauchy-Riemann operators. We start with the 
Dirac operator and consider the disturbed operator 

D + a = y^e 7 - h a 

^— ' OX-j 

j= i 

and 



(D + a)(D - a) = DD - a 2 = -A - a 2 , 
where a is a scalar-valued function. But we can do even more 

(D + a)(D + a) = DD + aa = —A — \a\ 2 . 
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An extensive study of these generalized or disturbed Dirac operators can be found 
in EDI and EH. 



2.5.5. Parabolic Dirac operators. Recently (cf. |H|), it was also possible to fac- 
torize the heat equation by a Witt-type basis. Which means to add two additional 
generating elements f and f + , where 

ff + + f + f = 1, 

f 2 = (f + ) 2 = o, 

fej + ejf = 0, 
f + e, + e,f+ = 0. 



and 



where 



j=i •? 



X,tJ 

i.e. D+ t factorizes the heat operator. The fundamental solution for the heat 
operator is 

, , H(t) _Ml 

where H(t) denotes the Heaviside-function. More on parabolic Dirac operators 
and the application to the non-stationary Navier-Stokes equation can be found in 
the already mentioned paper |8j. 

3. Parabolic Dirac operators related to the Schrodinger equation 

Let us denote with C the space-time domain C = G t x [0, oo) C {(x, t) e 
R n x R + } and by dC = G U M, where M = T t x [0, oo), its sufficiently smooth 
boundary. With the subindex t we emphasis the fact that the space domain (and 
therefore M) can depend on the time level t, but has to be bounded for all times 

t e [0, oo). 

The consideration of the Schrodinger equation is similar to the heat equation 
but the Schrodinger equations reads as 

-id t - A„, 

where A„ again denotes the Laplacian in M n . It is easily seen that we obtain a 
fundamental solution of the Schrodinger equation from the fundamental solution 
of the heat equation by the transform t — > it, hence 

V ; (2V^it) n 
where H(t) denotes the Heaviside-function. 
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We consider the operators 

j=i 

and hence 

{D% t fu = D% t D% t u = (D x + fdt ± if)(D x u + fdtu ± *f + «) 
= AcAeU - f£> x d t u =F if + D x u + f.D x d t w 
+ f 2 d 2 u ± ff + d t u ± if + D x ti + ±if + fd t u - (f + ) 2 u 

= -Au ± z(ff + + zf + f)<9 t w = -Am ± iflfeu. 

Next, we consider the reduced differential operator 
If we introduce the sigma-form 

dcr X) t = D Xit \ dV x dt 
we write Stokes theorem as follows 

/ / da x>t g = [ (/ <? + / (D x , t g) dV x dt 

JdC JC 

and thus 

/ / da x>t g = [if D+ t ) g + / (D~ t g) dV x dt. 

JdC JC 

Prom the fundamental solution of the Schrodinger equation we can switch to the 
fundamental solution of —A + id t by making a reflection t — > — t and applying D x t 
from the right we obtain the fundamental solution 

e(x, t) = E(x, -t)D+ tt = - P=- (D x + fd t + if + ) 

(2^/m(-t)) n 



H(-t) ,M 2 -i ^ c n ix 2 \ ..I 

/^FiH))« y 2 t^ J] ' V 2t At 2 J 1 



which gives the Borel-Pompeiu formula for the operator D~ t : 

/ e(x-xo,t-t )da Xjt g(x,t) = g(x Q ,to)+ e(x - x , t - t )(D~ t g) dV x dt 
JdC Jc 

and for g e ker D~ t Cauchy's integral formula: 

e(x -x ,t- t ) da X)t g(x,t) = g(x ,t ). 

ac 
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4. FACTORIZATION OP THE SCHRODINGER EQUATION 

The factorization of the Schrodinger equation is treated in several papers. Mostly 
the factorization is related to systems of differential equations in mathematical 
physics PJ, [Hj, [Hj. In [19j especially the relationship to the Riccati equation 
is investigated and quaternionic generalizations of the Riccati equation had are 
established. 

A detailed factorization of the stationary Schrodinger equation for n = 3,4 in an 
quaternionic context was treated in |3j. 

The Miura transform can be obtained by factorizing the ID Schrodinger equation: 

;U — v(x)u = — ( - — h a(x)^ ( a(x) N ] u. 



We have 



dx 2 \dx ) V dx 



— + a(x) ) ( — — a(x) ) u 



and thus 



dx J \dx 

( d 2 u d du 2 / \ 

[d^-Tx {a{x)u) + a{x) Tx- a{x)u . 

( d 2 u ( da o, ,\ 

= -{dx- 2 + {-dx-- aix) u 



- — ha (x) = —v(x) 
dx 



4.1. Factorization of the stationary Schrodinger equation. The Helmholtz 
equation can be treated quite similar to the Laplace equation and can be used to 
describe Maxwell's equations. The relationship between Maxwell's equations as 
a system of first order differential equations and the formulation in terms of the 
Helmholtz equation is nicely characterized by the factorization of the Helmholtz 
equation. The Helmholtz equation and the factorization in generalized Dirac op- 
erators is discussed in [20j, [21], [THj and [26j. In higher dimensions it is quite clear 
that the Helmholtz operator may be factorized by using Dirac operators: 

-A- k 2 = (D + k)(D-k). 

More difficult is the case of a variable potential V (x). We consider the case u = 
u(x) = u(x)e and look for suitable functions a with 

(-A - V {x))u =(D + a{x)){D - a(x))u 

= DDu — D(a(x)u) + a(x)Du — a 2 (x)u 

= —Au—Du ■ a(x) + a(x)Du — Da(x) ■ u — a 2 (x)u. 
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The underlined part does not vanish, because of the non-commutativity of the 
multiplications of elements of a Clifford algebra. Thus, we will change our ap- 
proach using a multiplication operator M a ^ defined by 

M a{x) u(x) := u{x)-a{x). 

Therefore 

(-A - V (x))u =(D + M a{x) ) (D - M a{x) ) u 

= DDu — D(ua(x)) + Du ■ a(x) — uDa(x) — ua 2 (x) 
= —Au — (Da(x) + a 2 (x)) u 

or 

Da(x) + a 2 (x) = V (x). (3) 

Equation (0) will be called Miura transform. It is nonlinear Clifford- valued first- 
order partial differential equation or equivalently a nonlinear real- valued first-order 
system of partial differential equations. 

4.2. Factorization of the non-stationary Schrodinger equation. Now, by 
using the operator D~ t we can also obtain a factorization of the non- stationary 
Schrodinger equation. Let us start with 

(D- t + M<^) {D- t -M<^)u(x,t) 

= D x,t D x,M x ^) - AmKMMM)) + ( D x,t u ( x ^)) a ( x ^) - u(x,t)a 2 (x,t). 
We have to evaluate 

D~ t (u(x, t)a(x, t)) = (D x + fd t - (u(x, t)a(x, t)) 
= D x (u(x, t)a(x, t)) + f(d t (u(x, t)a(x, £))) — if + (u(x, t)a(x, t)) 

If we again assume that u = u(x,t) = u(x,t)e is a scalar-valued function we can 
conclude 

= (D x u(x, t))a(x, t) + u(x, t)(D x a(x, t)) + f(d t u)a(x, t) + u(x, t)f(d t a(x, t)) 

—if + (u(x, t)a(x, t)) 

and finally 

(D~ t + M<^) {D~ t -M^)u(x,t) = 
= D~ t D~ t u(x, t) — (D x u(x, t))a(x, t) — u(x, t)(D x a(x, t)) — f(d t u(x, t))a(x, t) 
— u(x, t)f(d t a(x, t)) + if + (u(x, t)a(x, t)) + (D~ t u(x, t))a(x, t) — u(x, t)a 2 (x, t) 
= {—id t — A n )u(x, t) — u(x, t) (D x a(x, t) + f(d t a(x, t)) + a 2 (x, t)) 
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If we compare the last relation with the non-stationary Schrodinger equation we 
obtain 

D x a(x, t) + f(d t a(x, t)) + a 2 (x, t) = V(x, t) (4) 
This equation is similar to the equation 

d t g + Dg+\g\ 2 = (5) 

which has formal similarities to the canonical Riccati equation. Equation (0 was 
investigated in [16]. We specify the equation (J!J) in the following way. In appli- 
cations the potential usually depends only on x and not on t. Therefore, if we 
consider a potential V = V(x) equation (Jlj) suggest to consider also functions 
a = a(x) which depend only on x. Which implies that we will consider vector- 
valued functions a(x) depending only on x and we end up with the equation 

D x a{x) + g 2 (x) = V(x), (6) 

which is the same as for the factorization of the stationary Schrodinger equation. 

5. The nonlinear Dirac equation 

In this section we want to solve equation (JUj) under the following assumptions: 

(1) the potential V — V(x) is a scalar- valued function that depends only on 
x E G, where G is a compact and smooth domain of R n . 

(2) we are looking for a vector-valued function a(x) = Y^=i a i( x ) e i that de- 
pends only on x G G. 

If we apply the T-operator and use Borel-Pompeiu's formula (cf. Corollary EJ) 
we obtain 

a-Fa = T(V+ |a| 2 ), 
because we are only looking for a G im Q we have Fa = and we end up with 

a = T(V + \a\ 2 ) 

Our iteration procedure reads now as follows: 

a n := T(V + la^J 2 ), n=l,2, 

Obviously, V(x) + |a n _x| 2 G C reproduces a n = T(V + la,,^! 2 ) G G? ,n and 
because of tr T G imQ, also tr a Ln G imQ. 

That means that we also have to assume tra n-1 G imQ. 

Now, we consider the regularity. The Teodorescu transform is a weakly singular 
integral operator and therefore 

V g L P (G) ^TV G Wp(G), 1< p < oo, 
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on the other hand Sobolev's embedding theorems WHG) L P *(G), where p* < 
< n, with embedding constant C leads to the following. Set p* = 2p and 
assume that G W^(G). Then 

i/p 



\ VP 



-ll \\LP 




n n , „ v 



i/p 



3=1 3=1 

i.e. a n ^ G M^(G) implies K^l 2 G L P (G) and finally also T\a n _ x \ 2 G 
Thus we assume that n > p and 

np 22 n 

> 2p > 2np — p •<=>- 2p > np -<=>- 2p > n p > 



n — p ~ ~ 2 

Therefore we need the following assumptions too: 

(1) ra > p > | > 1, if n > 2 (or > p > 1 if n = 2), p G E, 

(2) a n _! G and tra n _ a GimQn Wp~*(T). 

The last assumption involves that if a n _ 1 belongs to WHG) due to the trace 
theorem for Sobolev spaces tra^^ belongs to W p P (T). 

In order to apply fixed-point theorems we first prove the boundedness of the 
sequence {a n }. We have 

HaJlvi/, 1 < \\T\\[ L P, (I 1^1 Up + II I«i„-i| 2 |Up) 
< ^ (ilV^lU, + C7 2 ||a h _ 1 ||^) 

<h\\V\\LP+k 2 \\a n ^\\ 2 Wpl , (7) 

where fci = ||T||[ iPi k 2 = hC 2 , and C is again the embedding constant for 

Wl ^ L 2 p. 

Lemma 1. We assume that ||V||lp < and denote by 



Then 
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implies that \\a +1 \\ w i < \\a \\ w i and the sequence {a n } n > no will be bounded 



by \\a nn \\wi- 



Proof. The inequality (JHD can be rewritten as 

-W < \\a„ \\wi - — <W 



-no 



p 2k 



2 



and hence is equivalent to 



and we obtain that 



2 1 k\ 

QnoWwl ~ T \\^n \\W£ + — \\V\\lp < 



and using (jZj) 

ll^no+lll^p 1 - ^lll^ll^ + fallOnJI^-i < ||o no ||vin- 



□ 



Theorem 1. Suppose that n > p > § > 1, pel, and \ \V\\lp < jj^jr- Then the 
equation 

a = T(V+\a\ 2 ) 

has at least one solution with 



l l h 

m\wi < 7TT- + 1 77 o - 1- \\v \\Lp- 



p 



2k 2 ' A/ Akl k 2 



Proof. Prom [T] we immediately see that 



K-i\\w}<± + W*\\a n \\ V Q<± + W. 



If we start with a e Wp(G) such that ||a ||wpi < 2^+^ then ^ ne sequence {a n } n eN 
is bounded from above by ^+W. Hence there exists a subsequence {a n ,} C W^G) 
with a n , -^aasn'^oo. Due to the continuity of T : L P (G) — > WHG) we conclude 

a = T(V + |a| 2 ). 

The norm estimates follows from the weak convergence of a n , in a convex set. □ 
But we can do better 

Theorem 2. We assume that 

(1) n >p > | > 1, pGM, 

(2) ||V||ii, < 

(3) a e WJ(G) and tra e imQ n W p 1_1/p (r), 

iiaiiwi < ^- - . /4t - iMiviiTT. 
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and define the sequence {a n } n eN by 

a n = T(V +\a n _ 1 \ 2 ), ra = l,2, .... 
Then there exists a unique solution a e WHG) with tr a G imQn wl~ l ^ p {T) of 

a = T(V+\a\ 2 ). 

Moreover, the sequence {a n }neN converges to a in and the solution a fulfils 
the estimate 



I II <r 1 1 fcl MT/ 



\LP- 



Proof. To apply Banach's fix-point theorem we investigate the contractivity of the 
mapping T(v + \a\ 2 ). At first we get 

WOin — %-lllWp 1 = ll^(|Oin-l| 2 - l«n-2| 2 )llwl < 1 1 |«„-l| 2 - I «n-2 1 2 I ■ 



litn-21 



Furthermore, 

1] \lp = \\(\a n -i\ ~ Ifi^jIXIfin-ll + K- 2 \)\\lp 

< II \<ki-l\ ~~ l^n-2l IU 2p ll ISi-ll + l^™-2 1 I U 2 ? 

^111 I _ I III III |2| |l/2 .mi 1 2 1 1 1/2 

^ II ISn-ll ISn-2l I \L 2 P\ ll&n-ll 1 1 LP > W \Qtn-2\ I 1 LP 

< II hn-l\ - K-2I IU 2p (C|l«n-lllwn + C'||a rl _ 2 ||iy p i 



^^Ilk-il-k^llUI^-J^-IlM 



LP 



Now, we use that a ra _ 1 and a n _ 2 are vector- valued and conclude 

111 11 111 ^ 1 1 i 111 ^ 1 1 i 1 2 1 1 1/2 ^ r^w 11 

II |«in-ll _ l«n-2l lli 2 f S II I%-1 — . ^n-2l \ \L 2 P S || l^-l" O n -2l || L P S U I l«n-l _ «n-2l I W£ ■ 

Then the contractivity constant can be bound from above by 



L = 2k 1 C 2 - W ) = 2k 2 - W) = 1 - 2k 2 W = 1 - Jl - 4fcifc 5 
V2/C2 / \2k 2 ) 



LP- 

□ 



Remark 3. Actually, it is no problem to find a suitable a e WHG) and 
tra e imQ n Wp~^ p {T). Obviously, a = fulfils all necessary conditions. But 
we have also other possibilities. Choose an arbitrary scalar-valued function 
b(x) such that \\b\\ LP is small enough then a = Tb is obviously vector-valued 
and fulfils all necessary conditions. 

Remark 4. More serious problem is that we had assumed the potential to be 
of small W-norm which need not to be fulfilled. In general the equation 

a = T(\a\ 2 ) + TV 
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is a system of nonlinear Fredholm equations of second with a weakly singular 
integral operator T. Hence methods to solve such types of nonlinear Fredholm 
equations could be applied. 

In what follows we will see that the solution of El is closely related to the solution 
of the linear Schrodinger equation. 

The following considerations are an easy generalization from the three dimen- 
sional case which was treated in [19] to the n-dimensional one. We go back to the 
Miura transform which could be considered as a Riccati equation as done in [19j. 
The Riccati equation 

du = pu 2 + qu + r, 

where p, q, r are functions, has received a lot of attention especially because of 
the wide range of applications in which it appears. This equation can be reduced 
into its canonical form 

dw + w 2 = —v. 

Also the Riccati equation is related to the one-dimensional Schrodinger equation 

- An - vu = (9) 
with a function v by the substitution 

du si 
w = — = omu. 

u 

Therefore the logarithmic derivative (in the sense of Marchenko) of a scalar- valued 
function u such that u ^ in G, is denned as 



du = u 1 Du = D(\n 



u . 



Proposition 1 ((cf. [19])). The scalar-valued function ip is a solution of the 
Schrodinger equation (G|) if and only if a := dip is a solution of the Miura 
transform 



Da = v + \a\ 2 . 



Proof. Suppose that there exists a function <p such that a = dip. Then 
Da = -^\Dp\ 2 - -Aip, and \a\ 2 = ^\D<p\ 2 , 

ip l tp> Lp Z 

and thus 



Da — \a\ 2 = Aw = v <^=^ —Aip — vtp — 0. 

Conversely, if (p is a solution of the Schrodinger equation a = dip is a solution of 
the Miura transform. □ 
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6. Application to the Gross-Pitaevskii equation 

Our considerations of the physical background of the Gross-Pitaevskii equation 
is based on the paper [31] and [32]. The Gross-Pitaevskii equation for the order 
parameter of the mean field theory is derived by the Fermi's zero-range pseudo- 
potential [lOj and [30]. 



6.1. Bose-Einstein condensate. An ideal Bose gas is a quantum-mechanical ver- 
sion of a classical ideal gas. It is composed of bosons, which have an integral value 
of spin, and obey Bose-Einstein statistics. The statistical mechanics of bosons 
were developed by Satyendra Nath Bose for photons, and extended to massive 
particles by Albert Einstein who realized that an ideal gas of bosons would form 
a condensate at a low enough temperature, unlike a classical ideal gas. This con- 
densate is known as a Bose-Einstein condensate. 



In dilute gases, only binary collisions are taken into account. The Hamiltonian of 
the interacting system is: 



H 



x 



dxdx'i>\x)i>\x')U(x - x')${x')${x), 



where \fr is the annihilation operator, m the mass of a single boson, H reduced 
Planck constant. Experimental realizations of Boson-Einstein condensate with 
dilute atomic gases has been achieved in a trapping potential V{x), whose shape 
is harmonic. 

The time evolution of the condensate is given by the Heisenberg equation 



iH^(x,t) 



tf, H 

2 



A + V(x) + / ^V, t)U(x - x')ty(x', t) 

2m J 



V(x,t). 



The mean-field approximation consists in replacing ^ by the classical field In 
the integral containing U(x' — x), this is a bad approximation for short distances. 
However, only collisions at low energy are relevant, which can be described by a 
single parameter: the s-wave scattering length a. Then one can write: 

U(x' — x) = g5(x — x) with g = 
And we obtain the Gross-Pitaevskii equation 



rn 



ih^ (x,t) 



H 2 

-—A + V{x) + g\y {x,t)\ 2 ) V (x,t), 



which is valid if 
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• the number of atoms in the condensate 3> 1, 

• low temperature 

• a <C average distance between atoms: vo? <C 1 {y is the density of the 
gas). 

Here, we have the Fermi pseudo-potential (given by the Dirac distribution) has 
one parameter - the scattering length. Thus, the simplest model which can de- 
scribe nonlocal repulsive interactions involves introduction of an effective potential 
with two parameters: the strength of interaction and the interaction range. The 
symmetry properties and the assumption of short-range interactions is enough to 
single out the effective potential with two parameters. For spherically symmetric 
interaction the U depends only on the relative distance U(x — x') = U(\x — x'\). 
Its Fourier image is then a real function of the squared wave number k : 

U = U(k 2 ). 

For the Fermi pseudo-potential is U = 1. In the next order of approximation we 
take the Lorentzian: 

U=(l + a 2 k 2 y 1 . 

The Lorentzian has been used recently as an approximation of the finite-range 
potential. In three dimensions the Lorentzian corresponds to the Yukawa effective 
potential: 

U 3D = -: \ exp (-- 

and in two dimensions the effective interaction potential is given by 

1 ( x 

where K (z) is the MacDonald's function. 

In the following we will use the operator form for the effective potential: 

A = 1 - a 2 V 2 = 1 - a 2 A. 

Which allows to consider the nonlocal Gross-Pitaevskii equation as a system with 
an additional dependent variable F describing the nonlinear term: 

(l-a 2 V 2 )F=h/f, (10) 

lh Ft + ^ A " v (x) " 9F ) ^ o(x ' t] = °' 

6.2. Stationary Gross-Pitaevskii equation. 

We obtain the stationary Gross Pitaevskii equation by setting 

V (x,t) = ip(x)e- i " t/n , where \i is the chemical potential, 
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and then 

h 2 \ 

~^A + V(x) + g\ip(x)\ 2 J ip(x) = fup(x), or 

.2w2 



l-a 2 V 2 )F(x) = \^(x)\ 2 , (12) 



2 



fv 

-—A + gF{x) + V(x) - fi \ ip{x) = 0. (13) 

Where the nonlinear equation is equivalent to a system of linear equations. Here, 
we see that the stationary Gross-Pitaevskii equation fits perfect in our schema. 
Equation (fT2j) is a Helmholtz equation 

with constant wave number and equation (JT3*j) is a (stationary) Schrodinger equa- 
tion 

—2m \ 
-A - — (gF(x) + V(x) - cp(x) = 0, 

of which ip is a solution if and only if 

a = dtp = (p~ 1 Dip = D(lmp) 
is a solution of the Miura transform 

— 2777, 

Da - \a\ 2 = -— (gF(x) + V(x) - (i) . 

6.3. Non-stationary Gross-Pitaevskii equation. 
Also the non-stationary Gross-Pitaevskii equation 

d ( h 2 

ih-ip(x,t) = i-—A + V{x)+g\i;{x,t)\ 2 

maybe written as a system of linear equations: 

(l-a 2 V 2 )F(x,t) = \i:(x,t)\ 2 , (14) 

~ lH §i ~ 2^ A ) + Vi - X) + 9Fi - X ' £) ) f) = °' (15) 



Equation IHM is again equivalent to a Helmholtz equation with constant wave 
number but the second equation (fT5|) has due to the part F(x, t) a potential which 
depends on the time t. 
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